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ABSTRACT 



< 

' Observations of the interstehar medium (ISM) and molecular clouds suggest these 

(-H ' astrophysical flows are strongly turbulent. The main observational evidence for tur- 

0H' bulence is the power-law energy spectrum for velocity fluctuations, E{k) oc /c", with 

O ■ £ [—1-5, —2.6]. The Kolmogorov scaling exponent, a = —5/3, is typical. At the same 

, time, the observed probability distribution function (PDF) of gas densities in both the 

Ci ', ISM as well as in molecular clouds is a log-normal distribution, which is similar to the 

initial mass function (IMF) that describes the distribution of stellar masses. 
^ ^ In this paper we examine the density and velocity structure of interstellar gas tra- 

' versed by curved shock waves in the kinematic limit. We demonstrate mathematically 

' that just a few passages of curved shock waves generically produces a log-normal density 

■ PDF. This explains the ubiquity of the log-normal PDF in many different numerical 

. simulations. We also show that subsequent interaction with a spherical blast wave 

, generates a power-law density distribution at high densities, qualitatively similar to 

I the Salpeter power-law for the IMF. Finally, we show that a focused shock produces a 

downstream flow with energy spectrum exponent a = —2. Subsequent shock passages 
reduce this slope, achieving a —5/3 after a few passages. We argue that subsequent 
^ • dissipation of energy piled up at the small scales will act to maintain the spectrum very 



near to the Kolomogorov value despite the action of further shocks that would tend to 
reduce it. These results suggest that fully-developed turbulence may not be required to 
explain the observed energy spectrum and density PDF. 

On the basis of these mathematical results, we argue that the self-similar spherical 
blast wave arising from expanding HII regions or stellar winds from massive stars may 
ultimately be responsible for creating the high mass, power-law, Salpeter-like tail on an 
otherwise a log-normal density PDF for gas in star forming regions. The IMF arises 
from the gravitational collapse of sufficiently overdense regions within this PDF. Thus, 
the composite nature of the IMF — a log-normal plus power-law distribution — is 
shown to be a natural consequence of shock interaction and feedback from the most 
massive stars that form in most regions of star formation in the galaxy. 
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Subject headings: ISM: structure, shock waves, stars: formation, stars: mass function, 
turbulence 



1. Introduction 



The distribution of stellar masses, known as the initial mass function (IMF), is of paramount 
importance to many fields of astrophysics. The form of the IMF plays a central role in subjects as 
diverse as galactic evolution and the formation and evolution of exoplanetary systems. Although 
its origin has been the focus of many analytical models for several decades, only recently have 
numerical simulations becom e available that cari include many of the impo rtant physical processes 



involved (see the reviews by iMcKee fc Ostriker 



2007 



Bonneh et al.l 120071 '). The form of the IMF 



is variously described as a log-norr nal distribution at low stellar masses with a po wer-law tail a t 
masses exceeding a solar mass (e.g. IChabrierll2003l ) . or as a multiple power-law (e.g. lKroupall2002l ). 
The Salpeter power-law index of -1.35 for the high mass power-law tail appears to be universal. 
These properties of the IMF must ultimately reflect robust properties of the dense substructure 
within molecular clouds in the interstellar medium (ISM), in which stars are born. 

It is not surprising therefore that one of the central problems in star formation is to characterize 
the mass distribution of the dense star- forming regions within clouds. Millimetre and submillimetre 
wave observations show that these clouds are highly inhomogeneous and are dominated by systems 
of filaments, punctuated by smaller denser regions in which clusters of stars form. Individual 
stars form in dense regions (n > lO^cm"^ ), whose mass distribution is known as the core mass 
function (CMF). The gas velocities in these clouds are observed to be supersonic and chaotic. 
Numerical studies reveal that "turbulent" supersonic gas motions can reproduce many aspects of 
this structure. A number of observational surve ys have shown that th e CMF can be modelled by a 
log-normal distribution in many instances (e.g. iGoodman et al.l 120081 ). Other studies suggest that 
the high mass tail of t his distribution is closer to a power-law, whose index is nearly identical to 
the Salpeter value (e.g. iMotte et al.lll998l : lJohnstone et al.ll2000l ). However, the de bate continues as 



to th e exact form of the mass distribution of these dense gaseous structures (e.g. iGoodman et al 
20081 ^. 



This subject has made major advances largely due to the advent of numerical simulations. 
These show that the filamentary nature and the mass spectra of structu re in clouds naturally resu lt 
from the action of supersonic turbulence within them (e.g., the review of lMac Low &: Klessenll2004l ) . 
The precise origin of the turbulence is still somewhat unclear. Molecular clouds are themselves 
embedded in a multi-component ISM and could be formed by several processes, including cloud 
collisions, spiral shock waves, and a combination of gravitational and magnetic instabilities. They 
are also shocked by the effects of feedback from star formation itself including expanding HII regions, 
supernovae, and powerful winds from massive star s. Simulations have started to make substantia l 
progress in following most of these processes (e.g. 



Wada &: Norman 



2001 



Tasker fc BrvanI l200ffl . 



The presence of supersonic turbulence within all molecular clouds has been interpreted as evidence 
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that they are short hved structures that are dissipated in one or two shock crossing times: star 
formation oc curs in the few transient turbulent stru ctures that have sufficiently high density to 
collapse (e.g. lElmegreenl l2002l : iHartmann et al.ll200ll ). Indeed, observational studies of molecular 



clouds in the nearby LMC galaxy indicate that the time scale to form star clusters is rapid (of 
order a Myr) an d that molecular clouds are rapidly dissipated in a few Myr as a consequence of 
their formation (jFukui et al.lll999l ). Turbulence within such clouds may result from processes on 
larger scales (> lOOpc) that tap into energy released by galactic shear, gravitational instability or 
large-scale expanding HII regions. 

In spite of the wide range of conditions in the interstellar medium and the details of the mod- 
elled process, time-dependent numerical simulations show that the density structure can be well 
modelled by a log-no rmal PDF over several orders of magnitude: from the diffuse atomic gas to 
molecular clouds (eg. IWada &: NormanI 120071 : iTasker BryanI 120081 ). At later tim es in their simu- 
l ation s of the ISM undergoing feedback from the effects of massive star formation iTasker &: Bryan 
(|2008l ) found that a power-law fit might also be possible. Numerical simulations of density fluctua- 
tions in purely isothermal supersonic turbule nce have a log-norma l PDF , and this is often taken as 
evidence for turbulence in the ISM. However, iPadoan Nordlundl (j2002l ) used the central limit the- 
orem to show that a flow with a power-law energy spectrum will necessarily have a log-normal PDF 
of density. Thus, a log- normal PDF of density does not provide any more evidence for turbulence 
than a power-law energy spectrum. 

How does the distribution of stellar masses arise from this density structure that characterizes 
the ISM? Stars are formed in fluctuations whose density exceeds a threshold for gravitational 
collapse. Recent theoretical work has emphasized that the shape of the IMF may be a combination 
of a power-law at large mass scales, that transitions to a log- normal form at lower masses. The peak 
of th i s distribution is a niass c haracteristic of gravitational collapse (e.g. iHennebelle &: Chabrier 



2008 



Padoan Nordlundll2002l ). The power-law tail and near Salpeter index for the mass function 



in molecular clo uds has been modelled as arising from shocks generated in nearly Kolomogorov 
turbulence (e.g. iPadoan k. Nordlundl l2002l ) . where the power-law index depends on the spectral 
index of the turbulent flow, and the half-width of the log-normal distribution depends on the 
turbulence Mach number. A recent analytic approach argues that the log-normal distribution for 
the IMF at low masses arises from that part of the gas density distribution that is supported 
by thermal pressur e, whereas the Salpeter tail a rises for higher mass cores that are supported by 
turbulent pressure (jHennebelle Chabrieil[200^ ') . An alternative explanation for a joint log-norma[ 
plus power-law distribution initial distribution is that an initial log-n ormal distribution ca n develop 
a power la.w tail if cores accrete over a distribution of time scales (jBasu &: Joned |2004| ) . Finally, 
ElmegreenI (|2002l ) noted that if one assumes that all gas above a density threshold in a log-normal 
distribution can form stars, then it is possible to r ecover the well-known S chmidt law that governs 
the global star formation rate in galaxies (see also lWada &: Normanll2007l ). . 



In this paper, we examine the mathematical properties of shock-driven gas motions and propose 
a new approach to explain the nature of gas motions and density structure in the ISM and molecular 
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clouds. We use analytical theory to examine both the density distribution expected in the gas due 
to the interaction of curved shocks, as well as the nature of the velocity fluctuations downstream 
of a shock. We show that the passage of just a few shock waves can very quickly establish log- 
normal density distributions. The passage of a spherical shock (i.e. blast wave) adds a power-law 
tail at large mass densities. We demonstrate that in general a log-normal distribution with a 
power-law at large densities is expected in media that are occasionally traversed by such large scale 
spherical shocks. These spherical shocks have been observed in the expanding HII regions that are 
a consequence of massive star form ation. Thus, although gravity may be important for the large 
scale dynamics of molecular clouds (jGoodman et al.l 120081 ) . the large scale shock waves of varying 
symmetry play an essential role in shaping the mass function of the cores. 

The motions induced in the wake of curved shocks are vortical in nature. Power is distributed 
to vortical motions across a wide range of scales without a cascade process that is essential for 
Kolomogorov turbulence. Although tempting, we will argue that it is problematic to interpret the 
observations of Kolmogorov-like scaling in terms of hydrodynamic (or MHD) turbulence. 

The paper is organized as follows. In ^we contr a st the properties of shock-driven and Kol- 



mogorov turbulence. Then in ^ we review iKevlahanI (jl997l )'s theory for vorticity generation by 



shocks propagating in nonuniform flows. Although it is co mmonly thought t hat straight shocks, 
weak shocks, and spherical shocks do not generate vorticity IKevlahanI ()1997l ) demonstrated that 
this is not true for shocks in inhomogeneous flows. In addition, we highlig ht the fact that curved 
shocks eventually focus. This focusing produces a pair of shock-shocks (jWhithaml Il974l ) which 
generate vortex sheets downstream of the shock. This effect appears not to have been considered 
before in astrophysical flows. In fjHwe show how multiple shock interactions could generate the ob- 
served log-normal and power-law distributions of mass density, and in f|5]we show that the observed 
velocity energy spectra could be produced by the quasi-singular vorticity generated downstream of 
focused shocks and blast waves. Finally, we interpret the results in terms of an astrophysical model 
for the role of shock waves in generating density structure in the ISM and molecular clouds, and 
its connection to the IMF and close with our conclusions for shocks and star formation in the ISM 



2. Shock-driven flow and Kolmogorov turbulence 



Observations of density and velocity fluctuations have suggested that many astrophysical flows 
are strongly turbulent. This phenomenon is widespread and includes a diverse set of systems — 
including H I emission in the interstellar mediu m (ISM), interstellar scin tillations (ISS), 100/i 
IRAS emission in the interstellar medium (ISM) ([Elmegreen &: Scald 12004 ) . vel ocity and density 



struc ture in molecular clouds, H I emiss ion in the large Magellanic cloud (LMC) (jElmegreen et al 



fluctuations in the solar wi nd (jHorburvl Il999l : iNicol et al.l 120081 ) . and hot II2O emission 
in accretion disks (jCarr et al.l 12004 ). The evidence suggesting that these fluctuations are in fact 
turbulent is principally their turbulent-like velocity dispersions, or, equivalently, energy spectra and 
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second-order structure functions, as well as their spatial density structure. S ometimes the scaliri g 
of higher-order structure functions is also taken as evidence of turbulence (e.g. iPadoan et al.ll2003l ). 
However, these results are not conclusive since there is still no rigorous th eory for the scaling o f 
high-order structure functions (only models, such as the one proposed bv IShe k. Levequd Il994l ). 
and there is not enough data to obtain properly converged statistics for structure functions of order 
greater than about 6 or 7 without using the extended self-similarity (ESS) correction. 

One of the central ideas in Kolmogorov's classical theory of turbulence is that energy is in- 
jected at large scales and cascades without loss through intermedi ates scales, until it reach es the 
smallest scales where it is finally lost to molecular dissipation (see lElmegreen &: Scald |2004| . for a 
detailed review). Fully developed incompressible turbulence is characterized by a power-law energy 
spectrum E{k) oc /c" where a = —5/3 for incompressible homogeneous isotropic three-dimensional 
turbulence and a S [—3, —4] for the enstrophy cascade in two-dimensional turbulence. Astro- 
phys ical fluctuations have observed power-law spectra with a G [—1.5,-2.6] (lElmegreen Scale 
20041 ). suggestive of turbulence. Radio propagation observations of the diffuse ISM have found 
that density fluctuations obey the "Big power law in the sky": Kolmogorov-li ke scaling of th e 
energy s pectrum extends over 1 1 ord ers of magnitude, from 10'' cm to 10^^ cm (jSpangleii Il999l )! 
However, iLazarian &: PogosvanI (|2000l ) comment that "the explanation of the spectrum as due to a 
Kolmogorov-type cascade f aces substantial difficul t ies." Indeed, th ey emphasize that "The existence 



1995 



Spangleijll999l ) is one of the great astrophysical 



of the Big Power Law (see [Armstrong et al. 
mysteries." 

Measurements of the velocity dispersion of gas within molecular clouds are made as a function 
of the size of the region in which the dispersion is measured. The data sho w that there is a power-law 



relation between the linewidth and size of a region such that Aw oc ( Ballesteros-Paredes et al 



20071') . This scaling is a direct consequence of the energy spectrum since /3 = —(a -|- l)/2. iLarson 
(Il98ll ) first deduced that (3 ~ 0.38, which is close to Kolmogorov turbulence (/3 = 1.3, a = —5/3), 
and was the first to suggest that turbulence must play a very important role in the gas, and as a 
consequence, in the process of star formation that occurs i n such clouds. Recent surveys for whole 
giant molecular clouds (GMCs) find [3 ~ 0.5 — 0.6 (e.g . iHever &: Brunt l2004l ). while studies of 



clouds with regions of low surface brightness find (3 ~ 0.4 (iFalgarone et al 



19921 ). Compressible gas 



motions which characterize the ISM are damped very quickly in shocks 
time of the driving scale. 



typically in one crossing 



Turbulence can also be characterized by the scaling of Cpi the exponent of the p—th. order 
structure function. 



{\u{x + r) — u{x)Y 



oc r 



(1) 



Kolmogorov's theory (jPrischl Il995l ) predicts that Cp is a linear function of p (with = p/3). 
However, experiments show that (^p is in fact a concave function of p, increasing more slowly than 
linear with order p. Some atter npts have been in ade to measure structure function exponents for 
astrophysical fluctuations (e.g. iNicol et al.l 120081 ) . although lack of data restricts the analysis to 
relatively low order (p < 4 or 6) and quantitative comparison with turbulent flows is difficult since 
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there is no accepted theory for how should scale. 

The equation of spectral energy balance for approximately incompressible homogeneous de- 
caying turbulence is 

dtE{k) = -2uk^E{k) +T{k), 

where k = \k\ is the magnitude of the wavenumber, the first term on the right hand side is viscous 
dissipation of energy (active only at small scales) and T{k) measures the rate of energy transfer 
from all other wavenumbers k' to k due to nonlinear interactions, i.e. it quantifies the energy 
cascade. Thus, in order to demonstrate conclusively the existence of an energy cascade one must 
estimate the energy transfer function, 

T{k) = [ u*{k) ■ P{k){ux~lj{k)) dS{k), 

Jk=\k\ 

where ( ) denotes the Fourier transform, u is the vorticity, P{k) is the divergence-free projection 
(for approximately incompressible flow) and the integral is over spherical shells in Fourier space. 
However, it is impossible to estimate T(k) from the available the observational data since calculating 
u X uj requires pointwise measurements of all components of the velocity and vorticity. 

The interpretation of power-law scaling of the energy spectrum in terms of fully developed 
Kolmogorov turbulence in the ISM is problematic for several reasons. 



1. 



2. 



Kolmogorov scaling is associated with inco nipressible neutral flow, whe reas the ISM is believed 
to be strongly compressible and magnetic. ISridhar &: GoldreichI (119941 ) proposed a theory for 
anisotropic incompressible magnetohydro dynamic (MHD) turbulence which gives an energy 
spectrum kj_^^^ in directions perpendicular to the mean magnetic field. However, the theory 
is not rigorous and it does not apply to the solar wind. Other weak turbulence calculations 



find kj_ or k^ 
the range from 



-3/2 



and stationary constant flux solutions may have e xponents anywhere in 
1 to —3 depending on the asymmetry of the forcing ([Galtier et al.ll2002l ). 



The "Big power law in the sky" extends over a range of scales that include many different 
physical processes, including scales where the gas dynamics approximations of fluid turbulence 
are not valid. How can the same scaling be maintained across scales with very different 
physics? 

It is not clear where the energy sustaining the ISM turbulence comes from. Candidates 
include massive stellar winds, supernovae, expanding HII regions, galactic rotation via spiral 
shocks, sonic reflection of shock waves hitting clouds, cosmic ray streaming, field star motions, 
Kelvin-Helmholtz and ot her fluid instabilities, ther mal instabilities, gravitational instabilities, 
and galaxy interactions (lElmegreen &: Scald l2004l ). 



4. 



As pointed out bv lLazarian &: PogosvanI (|200d ). the damping rate of MHD turbulence is much 
faster than previously thought: about one eddy turnover time (as for neutral fluids). This 
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implies very large energy injection scales and efficient and frequent forcing in order to sustain 
the turbulence (since the energy cascade takes about one eddy turnover time) . It is useful to 
recall that hydrodynamic turbulence typically has constant or frequent forcing, e.g. vorticity 
generation via the no-slip boundary condition, or the mixing layer instability in jets. Vorticity 
generation occurs on a huge range of scales and does not require the lengthy process of an 
energy cascade. 



In numerical simulation s, on the other hand, turbulence is often generated via spiral shocks o r 



super novae explosions (e.g. iJoung Mac Low 



2006 



Piontek &: Ostriker 



2005 



Wada et al.l I2OO2I ') . 



This turbulence is necessarily limited to low Reynolds numbers, and the power-law scaling of the 
energy spectrum is present over a small range of scales (about a decade). Supersonic turbulence 
from these and othe r simulations of th e ISM have a scaling close to t he —2 associated with the 
shock discontinuity (IKritsuk et al.l 120071 : IVazquez-Semadeni et al.lll997l ). 



Mathematically, the energy spectrum of a field is determined by its strongest singularity. If 
a function /(x) has a discontinuity in the p — 1 order derivative (where p is an integer), then the 
energy spectrum of f{x) has the form of a power-law 



E{k) ~ k 



'2p 



For example, a field containing shocks (discontinuities in the velocity field) has E{k) ~ k 
singularity 'worse' than a discontinuity would be required to generate a slope shallower than - 



2. A 

-5/3 



for ex ample an accumulation of discontinuities around a given point, or a fractal (see iHunt et al 
199d ). Note, however, that the converse is not true: smooth fields can have a power-law energy 



spectrum (e.g. adding together Gaussian functions of just the right size and amplitude could 
produce a k~^^^ spectrum). 

The fact that a power-law scaling is observed over such a wide range of scales suggests that 
a singularity may be res ponsible, such as the shocks that are ubiquitous in astrophysical flows. 
Kornreich &: Scald ([2000) have proposed galactic shocks propagating through interstellar density 
fluctuations as a way of forcing (or "pumping") supersonic turbulence. We go one step further: 
shock-generated vorticity alone may be enough to explain the observed energy spectra. Dynamical 
turbulence is not required. 

As mentioned above, simulations of supersonic turbulence typically produce a —2 spectrum. 
This is not su rprising since flows wi th turbulence Mach number Mt > 0.3 will spontaneously gener- 
ate shocklets (jKida Orszagjll990l ). Since the Kolmogorov —5/3 spectrum is l ikely not associated 
with a singularity, the —2 scaling dominates. However, as pointed out by lElmegreen &: Scalo 
(j2004l ). the observed spectra are often shallow than —2 and "the proposed shocks themselves have 
not been observed in real clouds." In this paper we suggest a shock-based explanation of the ob- 
served spectra that addresses both these issues. First, we show that focused shocks generate vortex 
sheets, which means that the shock-driven flow may have a spectrum E{k) ~ k~'^ even when the 
shock is no longer present. Secondly, we show that multiple passages of curved shocks produces 
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velocity fluctuations with a spectrum shallower than —2. Focused shocks are not the only shock 
structures that produce power-law scaling: we also find that multiple passages of strong spherical 
shocks produce similar results. 



Dobbs &: Bonnelll (|2007l ) have recently proposed a similar shock-based explanation of the ve- 
locity dispersion (i.e. energy spectrum) in molecular clouds. They use full smoothed particle 
hydrodynamics (SPH) simulations of the gas dynamics equations to show that shocks propagating 
through non-uniform fractal gas generate velocity dispersions close to the observed scaling. Our 
study is complementary: we use the full analytic expression for the vorticity produced by a curved 
shock in nonuniform flow to show that multiple shock passages could produce the observed velocity 
dispersion (and PDF of density fluctuations). We identify baroclinic vorticity generation as the key 
term. In addition, we find that fractal (or even non-uniform) initial conditions are not necessary. 



3. Vorticity generation by curved and oblique shocks 



In this section we review iKevlahanI (jl997l )'s theory for the vorticity jump across a shock in 
non-uniform flow. We emphasize two effects that are usually neglected in discussions of vorticity 
generation by shocks: the creation of vortex sheets downstream of highly curved shock regions (e.g. 
shock-shocks), and the role played by non- uniformities in the flow ahead of the shock (which can 
lead to significant vorticity generation even by straight shocks). 

The general expression for the vorticity jump in the binormal direc tion b across an unsteady 
three-dimensional shock moving into a non- uniform flow was derived by iHayesI (|l957l ). 



5ujh = n X 



djpCr 
dS 



S{p-'){pCr)-\Ds\Js + CrDsn)5ip) 



(2) 



where n is the shock-normal direction, s is the tangential direction, and the binormal direction 
is given hy b = s x n. Note that both s and b are tangential to the shock surface, and the 
normal component of the vorticity is continuous across a shock. d/dS is the tangential part of the 
directional derivative, = C — ^ is the shock speed relative to the normal component of the flow 
ahead of the shock A, and Ds is the tangential part of the total time derivative. A similar expression 
may be derived for the vorticity jump in the tangential direction s. However, the expression usually 
taken for the vorticity jump is 



6lj h 



n X U5 • K 



dCr 

as 



(3) 



where p is the normalized density jump across the shock (the shock strength), K is the curvature 
of the shock and U5 is the velocity tangential to the shock in the reference frame of the shock. 
This expression was derived by Hayes from ([2]) by assuming that the flow ahead of the shock is 
uniform. 
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KevlahanI (|l997l ) re-derived the vorticity jump equation, taking full account of terms due to 
the non-uniform upstream, finding that the vorticity jump in the binormal direction b 



5u ■ h 



A* 



dCr 



l + ^i dS 



P_ 



Du 



+ 



C}. 1 



Dt 1 + fJ.p 



Vp 



s] + ^LJ -b, 



(4) 



together with a similar expression for the vorticity jump in the tangential direction s. If the 
upstream flow is isentropic then a^/pVp ~ —Du/Dt, and if, in addition, the upstream fiow is 
quasi-steady and we normalize by the stagnation sound speed oq then ^ becomes 



6u: ■ b 



P 



1 



1 + p dS Ms \l + p 



P- 



1 



dS 



+ U X u ■ s 



+ pu -b, 



(5) 



where Mt is the turbulent Mach number of the upstream flow. We will use ([5]) in the remainder 
of the paper. Equation ([5]) may be simplified further for strong shocks by using the approximation 
/i« 2/(7-1). 



The density jump is given by 



5p = pp 



M2 



1 + 1/2(7 - 



P, 



(6) 



where p is often referred to as the shock strength. 



The first term on the right hand side of ([5]) represents vorticity generation due to the variation 
of the shock speed Mg along the shock; by symmetry it is exactly zero for spherical and cylindrical 
shocks. Because shocks are nonlinear waves (unlike acoustic waves), Mg is larger in regions of con- 
cave curvature and smaller in regions of convex curvature (with respect to the propagation direction 
of the shock). This difference in shock strength increases over time and eventually causes curved 
shocks to focus at regions of minimum c urvature, develo ping a flat shock disk bounded by regions of 
very high curvature (often called kinks) (lKevlahanlll996l ). In lab oratory experiments shock foc using 
is obtained by reflecting a straight shock off a curved surface (see lSturtevant Kulkarnvlll976l . for a 
detailed discussion and pictures of shock focusing experiments). The discontinuous shock strength 
at the kinks is called a shock-shock. In the ISM shock focusing could arise due to reflection off 
density gradients (e.g. vertically stratifled structure in disks), or due to small variations in shock 
curvature in blast waves. 

As mentioned above, since the flrst term on the right hand side of ([5]) is approximately singular 
at the location of a kink, extremely strong jet-like vortex sheets develop in the flow downstream 
of the kinks (see flgure[T]). These vortex sheets themselves have an energy spectrum E{k) ~ k~'^, 
and generate turbulence exponentially fast via the Kelvin-Helmholtz instability. Note that this 
scenario produces a —2 spectrum in the downstream flow, even when the shock is no longer present 
(resolving the first objection mentioned in the introduction). In other words, the —2 spectrum is 
associated with the downstream fiow, not with the shocks themselves (as has been assumed in the 
past). 
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The second term on the right hand side of ([5|) is barochnic generation of vorticity due to 
the misahgnment of pressure and density gradients as the flow passes through the shock. This is 
the dominant term for vorticity production across straight or weakly curved shocks. It is also an 
important term in the case of multiple shock passages as it nonlinearly mixes the flow, redistributing 
energy amongst different length scales (similar the the quadratic nonlinearity of the Navier-Stokes 
equations). Note that vorticity may be generated across a straight shock even if the upstream flow 
is irrotational. 

Finally, the third term on the right hand side of ([5]) is the additional angular momentum 
generated by compression of the flow in the direction normal to the front (i.e. conservation of 
angular momentum). This terms simply moves the entire energy spectrum up by the factor fi 
without changing its form. 

The following sections use simple examples to show how multiple shock passages can generate 
density PDFs and energy spectra similar to what is seen in molecular clouds. We consider three 
generic shock types: weak eddy shocklets (which form spontaneously in supersonic turbulence), 
focused shocks, and strong spherical shocks (which model the blast waves generated by super novae 
explosions). 



4. The distribution of mass density 

In this section, we derive the density PDF of interstellar gas that results by the passage of 
various shock waves. We first demonstrate that a log-normal distribution is very rapidly established 
in a medium that is repeatedly lashed by multiple shock passages (§ 4.1). However, in § 4.2 we show 
that a power-law behaviour for a density PDF is expected for the passage of a perfectly spherical 
blast wave. Interstellar gas can therefore be regarded as being characterized by a log-normal density 
PDF, which from time to time develops a power-law tail at high densities due to the passage of a 
spherical blast wave from a nearby supernova, or an ongoing stellar wind bubble. This situation 
typifies the gas dynamics in regions of star formation. 



4.1. Rapid generation of the log- normal density distribution 



We first present a very simple explanation for the origin of the log-normal distribution of 
density commonly observed in isothermal turbulent flows. It is well-known that flows with Mt > 
0.3 spontaneously g e nerat e small, relatively weak and highly curved shocks, called 'eddy shock- 
lets' iKida &: Orszagj (|l99d ). It is therefore reasonable to assume that a region of space will be hit 
several times by shocklets of varying strengths. If we assume that the density is approximately sta- 
tionary between shock passages (i.e. that the density changes primarily due to shock compression). 
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then from equation ^ the density after n shock passages is 

n 

Wr-r^ = TTn -u ,,0') 



'(X 



l[il + ^,(^){x)) (7) 



where we have normahzed density in units of the initial uniform density pQ. Let us consider the 
shock strengths fi^^\x) to be identicahy distributed random variables. Then, since 
we can take the logarithm of both sides and apply the central limit theorem to the resulting sum. 
This shows that the logarithm of density is normally distributed, i.e. the density PDF follows a 
log-normal distribution, 



^(^)=7iw''H — — )■ 

Note that application of the central limit theorem to derive ([8]) requires only that the random 
variables log(l + p^^\x)) have finite mean and variance. 

One might think that it would take hundreds of shock interactions to converge to this log- 
normal distribution. However, if the PDF of density jumps is symmetric, then the rate of con- 
vergence is quite fast, 0(n~^/^). In fact, if the PDF of density jumps is uniform (a reasonable 
assumption) then, as shown in figure EJ as few as three or four shock passages generates a very 
good approximation to the log- normal distribution. 

In particular, if the PDF of density jumps fi is uniformly distributed in (0,2/(7 — 1)] (i.e. 
between its minimum and maximum possible values), then the logarithmic mean and variance of 
the log-normal distribution after n shock passages are respectively 

bi^ = 1ln2±l (9) 

9 n , 7 + 1 , ^ 

" = 12 <i°> 

There are two important aspects to these results. First, we see that both the mean and variance 
of the log-normal are largest for nearly isothermal gases (i.e. gases with 7 ~ 1). Molecular gas 
can be well described as isothermal up to densities of n ~ 10^ cm~^ due to efficient CO molecular 
as well as dust cooling. Second, the result shows that the logarithmic mean increases proportional 
to the number of shocks n — and that the width of the distribution measured by the variance 
also grows proportional to n. These tre nds have been observed in the density PDFs of simulations 



of molecular clouds (e.g., the review of iMac Low &: Klessenll2004l ). The overall amplitude of the 



distribution decreases with growing n simply because the integral of P{p) must equal unity for any 
n. These trends are shown in figure O Note in particular how few shocks are required to establish 
a log-normal distribution to high accuracy. In a self-gravitating medium such as molecular gas, 
eventually gravity takes over and the dense cores begin to collapse. 

The explanation for the log-normal dis tribution of density propose d here is even simpler and 



more general than that given previously by lNordlund &: PadoanI (jl999l ). In fact, our results could 
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Fig. 1. — Focused shock model: — shock strength, shock profile. Note the shock-shocks 

(discontinuities in shock strength) which generate vortex sheets in the downstream flow. 
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Fig. 2. — Convergence to a log-normal PDF of density after n = 2,4,6 shock passages. , 

log-normal distribution, — PDF of density. 
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explain iNordlund &: Padoanl (|l999l ) 's observation that in a numerical simulation of isothermal su- 
personic turbulence ". . . the high-density wing of the Log- Normal is established very early - soon 
after the first shock interactions." It is precisely these first few shock interactions that generate 
the log-normal distribution. 



4.2. The power-law distribution at large densities 

We have shown that interacting weak shocklets typical of supersonic turbulence quickly gen- 
erate a log-normal PDF of mass density. However, observations show that the PDF of mass has 
a power-law tail at high masses with an exponent near the Salpeter index of -1.35. Although the 
mass and density PDFs are not identical, this observations suggests that the density PDF should 
also have a power tail at high densities. We show here that such a power-law tail may be explained 
by the interaction of the log-normal density fields of supersonic turbulence with a strong spherical 
shock (i.e. blast wave). 

We ad opt the solution f or strong spherical shocks with sustained energy injection E(t) oc 



derived by iDokuchaevI (|2002l ). This solution generalizes the Sedov-Taylor self-similar solution for 
a point blast explosion modelled by p = (i.e. an instant shock) to permanent energy injection 
modelled by p = 1 (i.e. an injection shock). The first case corresponds to the instantaneous 
addition of energy to the ISM (as in a supernova explosion), while the second case corresponds to 
the continuous injection of energy (as in a massive stellar wind). The instant shock corresponds to 
the classical Sedov-Taylor solution for super novae explosions. 

We assume that the PDF of finding a particular value of gas density pi is proportional to 
the the space-time volume where the density exceeds pi, 

rKpi) 

P{p > pi) oc / R^{t)dt (11) 
Jo 

where R{t) oc is the radius of the spherical shock at time t and t{pi) is the time at which 

the density behind the shock is equal to pi. Using the relation Ms{t) oc R{t)/t, equation ([6]) can 
be inverted to find t(pi) oc p5/(2{-3+p))_ 

Using the definition of the PDF, we find that the PDF of density due to the interaction of a 
homogeneous gas with a spherical blast wave has the form 

Pi^p) = _i / R^(t) dt OC p-(17+p)/(6-2p) (^2) 

dp Jo 

(where we have re-labelled pi as p). Therefore, the density PDF is a power-law with exponent 
— 17/6 ~ —2.8 for an instant shock and —9/2 = —4.5 for an injection shock. These slopes are 



We would like to thank an anonymous referee for suggesting this space-time derivation approach. 
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significantly steeper than the Salpeter value for the mass PDF of —1.35. However, the actual 
relation between the density PDF and the mass PDF depends on the precise assumptions made 
about the scaling of clumps (i.e. mass equals density times a lengthscale cubed). Thus, one should 
not necessarily expect the same index for both the density and mass PDFs (although the power-law 
form should be robust). 

Mathematically, the PDF of density resulting from the interaction of a spherical shock with a 
log-normally distributed density field is simply the convolution of the PDF (jl2p with the log-normal 
distribution of density ([8]). This produces a PDF which is log-normal for small densities, and has 
a power-law p^(i'^+P)/(^~2p) high densities up to a maximum density proportional to 1/(7 — 1). 
Figure [3] shows the resulting PDF, where the initial log-normal PDF is the result of four shocklet 
passages with maximum shock Mach number M = 1.5 in a nearly isothermal gas with 7 = 1.05. It 
is important to note that because of the efficiency of CO and dust cooling, molecular gas remains 
essentially isothermal to density of ~ 10^ cm~^ — hence our choice of an adiabatic index near to 
a value of unity. 

Over time, the continuous production of shocklets will cause the PDF to revert to log-normal 
form, at the slower rate 0{n~^) (since the PDF is not symmetric). Thus, the presence of a power- 
law in the density PDF suggests that the flow has interacted recently with a blast wave (such as a 
super nova explosion). 



5. Energy spectrum 

5.1. The multi-shock model 

In addition to explaining the log-normal and power-law distributions of density, multiple shock 
interactions could also explain the observed power-law energy spectra (or velocity dispersion). It 
is important to emphasize that we consider the energy spectrum of the flow downstream of the 
shock. The velocity discontinuity associated with the shock itself gives E{k) oc k~'^, which will 
determine the power-law exponent of the energy spectrum only if the downstream flow has an 
energy spectrum equal to or steeper than -2 (or unless there are no longer any shocks present). 
We will show that this is not the case in general for multiple shock passages. To the best of our 
knowledge, the special effects of focused shocks on the downstream flow have not been examined 
before in the astrophysical context. 



The relevance of the case of multiple shock passages has been established bv lKornreich &: Scalo 



([20001), who found that the average time between shock passages in the ISM is "small enough that 



the shock pump is capable of sustaining supersonic motions against readjustment and dissipation." 

We use a semi-analytic approach to calculate the vorticity generation due to single and multiple 
passages of curved shocks. The vorticity jump is calculated using equation and the velocity 
and the required gradients of upstream quantities are calculated using the fast Fourier transform 
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Fig. 3. — Generation of a power-law PDF at large densities in a nearly isothermal gas by a spherical 
blast wave interacting with a log-normally distributed density field. is the initial log-normal 

PDF and the slopes -17/6 and -9/2 eorrespond to an instant shock and an injection shock respec- 
tively. Note that the upper limit of the power-law range is proportional to 1/(7 — 1), and thus we 
expect the largest power-law ranges for nearly isothermal gases, i.e. those with 7 fs 1. 
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(FFT) on a computational domain with periodic boundary conditions. The computational grid is 
256^ in all cases. The initial flow is assumed to be irrotational. 

We make the following assumptions: 

1. Frozen vorticity: flow evolution is due to the shock alone. The flow is approximately steady 
between shock passages. This is similar to the rapid distortion app roximation for strained 
turbulence, and it linearizes the problem. iKornreich &: Scald (j200d ) make a similar frozen 
vorticity assumption to neglect flow evolution during the shock passage. We deliberately 
ignore the internal dynamics of the flow: energy re-distribution amongst scales is due to the 
shock. 

2. Strong shock: the shock does not change due to interaction with the flow. 

3. Steady shock: the shock's shape and strength distribution are fixed. 

4. Random shock: The direction and phase of the shock are chosen randomly for each passage, 
and the results are averaged over many independent realizations. This is the same assumption 
we made in deriving the density PDFs. 

This semi-analytic approach i s extremely efficient numerically, and is similar to the kinematic 
simulation method for turbulence (jFung et al.lll992l : lElliott Maidalll995l ) In kinematic simulation 
the energy spectrum is specified, but the complex phases are chosen randomly. Thus kinematic 
simulation is accurate for quantities that depend on second-order moments of the velocity field (e.g. 
particle dispersion and energy spectrum). However, as its name implies, kinematic simulation does 
not accurately represent the dynamics of a turbulent flow. In particular, the linearizing assumption 
that is at the heart of the model begins to break down once a significant amount of energy piles 
up at smaller scales. This will play an important role in defining the steepness of the energy 
distribution, as we shall see later in this section, and which will be discussed in the next. 



The shock profile and shock strength profile for each case are described separately below. 



5.2. Focused shocks and shock shocks 

We first consider flow driven by a focused shock. As mentioned in the introduction, a focused 
shock is characterized by a flattened shock disk bounded by two shock-shocks (or a shock-shock 
ring in the three-dimensional case). The shock strength is (approximately) discontinuous at the 
shock-shocks, which generate a vortex sheet (with spectrum A;"^) behind the shock. Vortex sheets 
are linearly unstable via the Kelvin-Helmholtz instability, and generate turbulence very efficiently 
and quickly. 

We model the shock profile (j) and shock speed Ms by 

<f>{x) = asin^(A;ix) (13) 
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Ms{x) = Afo + 1 + - 



vr 




(2j - 1) 



cos ((2j — l)2kix) sinc^ 



(2i - l)2h 



n 



(14) 



The expression for Ms is simply the Fourier series for a square wave, where we have used the 
Lanczos cr-factor (the sine term) to remove the Gibb's oscillations. The number of terms N in the 
series is taken equal to the number of Fourier modes used in the spectral method. We use a similar 
expression for a two-dimensional shock z = (f){x,y), with wavenumbers ki and k2 in the x and y 
directions. Although this shock is simply a model (i.e. it is not the solution of the nonlinear wave 
equation governing shock motion), it captures its main qualitative features. The shock profile and 
shock strength are shown in figure [H 

We consider two initial conditions: a uniform flow (i.e. constant density and zero velocity), 
and an irrotational flow with a Gaussian energy spectrum E{k) oc exp(— fc^). The latter flow models 
the final decay regime of a turbulent flow, when viscous diffusion dominates. 

Figure IHa) shows the energy spectrum of the downstream three-dimensional flow after one, 
two and three passages of a focused shock with ki = k2 = 1, Mq = 6 with zero velocity initial 
condition. Note that the initial k~'^ scaling of the energy spectrum (due to the velocity discontinuity 
associated with the vortex sheet downstream of the shock-shocks) becomes gradually shallower 
with each shock passage. This redistribution of energy to smaller scales is due to the quadratic 
baroclinic terms depending on the (inhomogeneous) upstream flow in the vorticity jump equation 
([5]). Although the effect is entirely kinematic, these quadratic terms redistribute energy amongst 
scales in a way analogous to the quadratic nonlinearity in the Navier-Stokes equations. 

After three shock passages the energy spectrum has a scaling similar to k~^/^. More shock 
passages would produce an even shallower power-law. Figure[Hb) shows that the form of the energy 
spectrum is relatively insensitive to the choice of initial condition. 



We now consider the case of perfectly spherical shocks described in §4.21 Due to symmetry, 
the shock strength Ms is constant along the shock and therefore the first term in ([5]) is identically 
zero. This means that vorticity production is due entirely to the baroclinic and angular momentum 
conservation terms. Since the initial flow is irrotational, only the baroclinic term is active for 
the first shock, however subsequent shocks generate vorticity both bar o clinically and via angular 
momentum conservation provided they are not coincident with the first shock. 

For strong spherical shocks the vorticity jump dH) reduces to 



5.3. Spherical shocks 



Su ■ b 



7 + 1 
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dS 



+ UJ X u ■ s 



(15) 



and thus, provided the upstream flow is non-uniform, the vorticity jump is proportional to Ms{r). 
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If the upstream flow is smooth and irrotational, the energy spectrum of the downstream flow 
due to a single shock passage is simply the conv olution of the Fourier transforms of the singular 



shock strength Ms{r) given bv iDokuchaevI (|2002l )'s blast wave solution, 



and the gradient of the turbulent kinetic energy (e.g. /c exp(— A;^)). This gives E{k) ~ k~'^ for an 
instant shock and E{k) ~ k~^'^^^ for an injection shock as A: — > oo. The energy spectrum for a 
generic self similar shock corresponding to (jl6p is 

i?(A:) ~ A:-(3+W{i+P/2). (17) 

Recall that these are the energy spectra of the flow downstream of the shock. 
The singularity at r = is removed by using the following regularization 

where the parameter is set slightly small than the grid size and Ms{0) is set to ensure that 
Ms = 1 at the edge of the computational domain. The upstream flow is assumed to be irrotational 
and to have a Gaussian energy spectrum E{k) oc exp(— A;^). 

Figure [5] shows the energy spectrum of the downstream flow after one, two and three passages 
of a spherical instant shock and a spherical injection shock. Note that although the injection 
shock generates a much steeper energy spectrum after one passage {k~^^^^ compared to k^^ for 
the instant shock), after three passages both spherical shock flows have a spectrum with a slope 
close to —5/3. As with the focused shocks, more shock passages would further decrease the slope 
as energy is redistributed to smaller scales by equation p^ . This result is a consequence of our 
kinematic treatment. As we shall argue next, nonlinear processes will cut in to limit this evolution 
of the energy spectra. 



6. Astrophysical implications 

While it is well kno wn that curved shocks are an effective way of forcing a turbulent flow 
(jKornreich &: Scaldl200Cl ). we have developed a simple kinematical model which demonstrates that 
shock interaction alone may can produce energy spectra of velocity fluctuations and mass density 
distributions consistent with the observations in the ISM. Fully developed turbulence is not nec- 
essary. Our result has several important implications for the density structure of gas in the ISM, 
and for the formation of stars within these structures. 
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Fig. 4. Energy spectra for multiple passages of focused shocks: slope decreases with each shock 
passage as energy is redistributed to smaller scales. Then spectra have been normalized so that 
£^(1) = 1. (a) Uniform initial flow, (b) Gaussian initial energy spectrum. 




Fig. 5. — Energy spectra for multiple passages of spherical shocks: the slope decreases with each 
shock passage as energy is redistributed to smaller scales. The initial energy spectrum is Gaussian 
and all spectra have been normalized so that E(l) = 1. (a) Instant shock [p = 0). (b) Injection 
shock {p = 1). 
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6.1. Energy spectra 

Simulations have confirmed that the structure observed in the interstellar medium and molec- 
ular clouds probably derives from shock-driven processes. Our results show that these processes 
do not drive turbulence in the classical sense of a systematic cascade in energy from large to small 
scales. The observed E{k) k'^/^ "big power-law in the sky" may actually be entirely shock- 
driven, and not a signature of fully-developed turbulence. This could explain why the power-law 
extends over such a huge range of length scales (spanning several different physical regimes, from 
the diffuse ISM to molecular gas), since the power-law of a shock-driven flow is due to a singu- 
larity, and so extends over all scales (until the viscous cut-off). Shock driven vortical motions are 
generic in the ISM, and we have shown that a log-normal ordering of gas structure develops rapidly 
over all scales of the gas. We demonstrated that this rapid appearance of a converged log-normal 
distribution is an expected property of shocked flows and the central limit theorem. 

Our results also show that shocks should be much more efficient forcing the flow than has 

been appreciated previously. This is because the shock immediately distributes energy down to the 
smallest scales, according to a spectrum very close to the Kolmogorov -5/3 profile. Relying on a 
classical turbulent energy cascade is not feasible, since the time needed to transfer energy from the 
largest to the smallest scales (roughly one large-scale eddy turnover time) is far too slow. 

How in this picture can a turbulence spectrum close to -5/3 be universal if repeated shocks 
can continue to make it shallower? As noted previously, continued steepening assumes a purely 
kinematic mechanism without the limitation of nonlinear or viscous effects. As we pointed out 
in the introduction, although radio propagation observations find a turbulence-like spectrum close 
to -5/3, the full set of astrophysical observations find spectra in the range [—1.5,-2.6]. So, our 
physical model needs to be able to produce a range of spectra, not just the single universal slope. 
Thus, we need to be able to suggest why -5/3 is the most likely spectrum, as well as accounting for 
other slopes. 

There is in fact a natural limit to how shallow the slope can be. The slope could never become 
shallower than -1, since this would imply infinite energy (assuming an arbitrarily small minimum 
length scale). With each shock passage the slope of the energy spectrum becomes shallower as 
the shock redistributes energy to smaller scales (much as the nonlinear term of the Navier-Stokes 
equations does). In the kinematic limit this process would continue until the slope approaches -1. 
However, in reality once a sufficient amount of energy accumulates at the smallest scales energy 
dissipation by viscosity becomes significant. At this point, the linearizing assumptions of the 
kinematic model break down. Energy dissipation limits the continued transfer of energy to smaller 
scales, as well as driving an energy cascade to small scales (as in decaying turbulence). The energy 
cascade necessarily produces an energy spectrum that converges quickly to -5/3, since we now have 
the conditions required for the existence of a universal inertial range: a sink of energy at small scales 
separated from a wide and continuous range of active scales. The fact that the initial condition is 
already close to -5/3 means that we expect this adjustment to happen exponentially fast (as in the 
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Kelvin-Helmholtz instability, where the initial spectrum is -2). Thus, the shock forcing effectively 
transfers energy to smaller and smaller scales until dissipation drives a turbulence cascade which 
fixes the slope of the energy spectrum at the universal value of -5/3. 

The essential point is therefore that shock-generated vorticity will very quickly establish a 
power-law spectrum that is close to the Kolomogorov value. We argue that at this point, energy 
dissipation will tend to keep it there. 



6.2. Feedback and the IMF 

Another important result of our analysis concerns the appearance of a power-law tail for 
initially log-normal density PDFs that interact with spherical shock waves. The spherical symmetry 
means that a power-law will be imposed on the original log-normal distribution. 

Consider the typical situation in a molecular cloud where a star cluster has started to form. 
We showed that the shock-driven motions that dominate an ISM with supersonic velocities rapidly 
produces a log-normal density PD F. It is well established that a ty pical star forms as a member 



of a star cluster (e.g. reviews bv iPudrita l2002l : iLada k. Ladal l2003l ) . The collapse of the dense 



gravitationally unstable regions results in the observed IMF. Before this process has terminated, 
however, a cluster is strongly impacted by the approximately spherical shocks associated with the 
most massive stars that have already formed nearby. Observations show that most embedded 
clusters are adja cent to HII regions that are excited by massive stars in nearby clusters (e.g. 



ElmegreenI l2002l ). Most clusters show evidence that their formation could, in fact, have been 
triggered by the powerful shock waves associated with the expansion of such nearby HII regions. 
We showed that the passage of these shocks alters the initial log-normal density PDF into one that 
has a power-law tail. This feedback from massive stars will therefore change the form of the IMF, 
most likely by producing a power-law tail. 

Subsequent spherical shock passages will further modify the index of the power-law tail. We 
note, however that at most one or two passages would be expected since cluster formation is 
typically completed in a million years, which is roughly the crossing time for such an event. Thus, 
we conclude that feedback from massive stars is likely to leave a signature on the form of the 
density PDF in the gas which will carry over into the IMF (from those fluctuations that undergo 
gravitational collapse). Our analysis of how triggering may effect the form of the IMF given an 
initial gas density PDF will appear in a future paper. 

An important caveat that we have not yet discussed is the possible role of magnetohydrody- 
namical (MHD) processes in shaping the density PDFs. Extensive sets of simulations have shown 
that magnetic fields with energy densities comparable to gravitational self-energy certainly can af- 
fect the form of the density PDF. However, the observations are best matched with fields strengths 
that are much less than this critical value. The simulations that best match the observations 
of dense magnetized cores involve "supercritical" field strengths, i.e. fields with significantly less 
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energy density than gravitational and which are therefore subject to gravitational collapse (e.g. 
Padoan fc Nordlundl[l999l : Ixillev k Pudritzj [20071 : ICrutcher et ahlboosl l Thus, MHD effects in the 
bulk of molecular gas are probably far less pronounced than envisaged in earlier models of star 
formation. 

Finally, and as technical aside, we suggest that it may be more useful numerically to generate a 
flow with the correct energy spectrum over a very wide range of lengthscales than to try to simulate 
the full nonlinear dyna mics of a turbulent flow over a very small range of length scales (as is done 
currently). For example. Elliott &: Maidal (|l995l )'s method is able to efficiently generate a Gaussian 
random field with a k~^^^ energy spectrum over 12 decades in two dimensions using only about 
47000 computational elements (compared to 10^^ computational elements for a conventional non- 
adaptive approach!). In addition to being computationally efficient, this approach may actually be 
a better model of the hydrodynamics of the ISM. 



7. Conclusions 

Our results have important implications for the origin and evolution of density fluctuations, and 
particularly the CMF in molecular clouds. Our combined log-normal and power-law distribution 
arises because there are two natural kinds of symmetry to shocks — planar and spherical — that 
combine in a natural way. The fact that the power-law energy spectrum and log-normal distribution 
of mass density are also observed in the diffuse ISM strongly suggest that shock-generated vortical 
motions play a profound role on all scales in the ISM, and even in star formation within molecular 
clouds. Our specific conclusions are summarized below. 

1. Supersonic turbulence with Mt > 0.3 spontaneously generates relatively weak and short-lived 
'eddy shocklets'. A few passages of these shocklets is sufficient to generate a log-normal 
distribution of mass density. Thus, a log-normal distribution of density should be typical of 
supersonic turbulent flows, and is established very rapidly (the passage of just a couple of 
shocks will suffice). This is contrary to the usual assumption that enormous numbers of shock 
passages would be required due to the slow convergence to a normal distribution. 

2. A spherical blast wave interacting with a log- normally distributed density field produces 
a power-law distribution of density at large densities, qualitatively similar to the observed 
Salpeter tail for the IMF. The power-law range increases like 1/(7 — 1), and thus will be 
largest for nearly isothermal gases (i.e. those for which 7 ~ 1). Over time the power-law 
gradually decays to a log-normal distribution due to the action of the continuously generated 
shocklets. Thus, the presence of a power-law in the density distribution implies that the flow 
has interacted fairly recently with a blast wave (e.g. supernova explosion). 

3. A single strong shock passage can generate a relatively steep power-law energy spectrum over 
all length scales (e.g. k~'^ for a focused shock or for a spherical blast wave) due to the 



-23- 



singular structure of the shock strength. In the kinematic hmit that we have investigated, 
subsequent shock passages increase the total energy and reduce the slope of the energy spec- 
trum as the quadratic nonlinear term representing baroclinic generation of vorticity by the 
shock redistributes energy to smaller scales. Three shock passages suffice to produce an en- 
ergy spectrum close to k^^^^. Note that molecular clouds could not support more than a few 
large scale events, such as expanding HII regions, without being destroyed. 

4. We argue that the onset of energy dissipation that is expected when energy piles up at the 
smaller scales acts to limit the energy spectrum generated by shocks to the Kolomogorov 
value. 

5. The energy spectrum we find is that of the downstream flow, not that associated with the 
velocity jump of the shock itself (which has a k~'^ spectrum), i.e. we are not measuring the 
spectrum of the shock itself. 

We close by noting that, to our knowledge, this is the first time vorticity generation and mass 
clumping by multiple hydrodynamic, curved shocks has been quantified analytically. 

The authors wish to thank Ethan Vishniac for helpful discussions on the observations and 
origins of astrophysical turbulence. We thank an anonymous referee for very useful comments on 
our manuscript. REP also thanks Eric Feigelson for interesting discussions on the general nature 
of log-normal plus power-law distributions in statistics. The research of both NK and REP was 
supported by NSERC Discovery Grants. 
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